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to O.lGeV 2 ) and x (down to 10 -6 ), are described. Conventional leading twist 
NLO perturbative QCD gives an excellent description of all the new data with Q 2 
above a GeV, leaving very little room for either higher twists or higher logarithms. 
We summarize the current status of NLO fits, and the determination of the gluon 
distribution and the strong coupling constant from structure function data. Finally 
we consider some of the theoretical questions issues by the new data. 



Until a few years ago our knowledge of structure functions and derived 
quantities such as parton distributions came almost entirely from fixed-target 
experiments, using muon and electron beams. Now it is being complemented 
and extended by results from the HERA ep collider, especially in the regions 
of low x and high Q 2 . At this workshop both the fixed target and the HERA 
experiments have presented new data on F2, extending the kinematical range 
and reaching a precision of more than a factor twoJpetter than the results from 
the 1993 run shown at last year's meeting in Paristl The ZEUS data now reach 
Q 2 values down to 0.1 GeV 2 , and x values as low as 10~ 6 . New measurements 
of the ratio of the longitudinal to transverse photon absorption cross section R 
were shown by NMC and CCFR. A first glimpse of Ft, at low x was presented 
by HI. There were also new data on the charm quark contribution to F2. 

Already at last year's meeting it was clear that HERA data at low x 
but with Q 2 above a few GeV 2 are well described by conventional leading 
twist NLO perturbative QCD evolution from a flat or valence-like boundary 
condition. The increased precision of the new data have considerably reinforced 
this conclusion, putting tight constraints on the size of novel effects such as 
BFKL logarithms (i.e. higher orders) and parton recombination (i.e. higher 
twist). The apparent absence of such contributions poses various interesting 
theoretical problems, some of which were discussed at this meeting. 

On the other hand the serendipitous success of NLO perturbative QCD at 
small x encourages us to use it to extract parton distributions in the proton, 
and in particular the gluon content of the proton at low x. At large x the 
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gluon can be constrained with the help of the newly available inclusive jet 
data from the Tevatron and an improved understanding of the older prompt 
photon data. Both of these recent developments lead to new information on 
the strong coupling. 

Here we only discuss in detail data on unpolarized nucleon structure func- 
tions: data on diffractive and polarized nucleon structure functions was pre- 
sented in the sessions of other working groups, and discussed in joint sessions. 
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Figure 1: The kinematic region covered by proton structure function data, including the 
new measurements shown at this workshop. 
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1 The Structure Function Data 



The differential cross section for neutral current deep inelastic scattering is 
related to the three structure functions F 2 ,R = F 2 /2xFi and F 3 according to 

d 2 a ep^eX 47rQ ,2 y 2 



dxdQ 2 Q 4 x 



(1 - y + 2{1 + R) ) F ^ T y(l - l/Wfo Q 2 ) 



(1) 

Here Q 2 is the four momentum transfer, x the Bjorken-x value and y the energy 
transfer from the electron to the proton in the proton rest frame. The parity 
violating structure function Fa only becomes significant in the region where 
Z° exchange or Z a — 7 interference dominates, i.e. at large Q 2 . The effect of 
R is important at large y, typically y <^ 0.4. 

The new kinematical region co|tered by the HERA and fixed target data 



is shown in fig|T[ The HERA dataHa were collected in 1994 and represent a 
tenfold increase in statistics with respect to the 1993 data, allowing a significant 
extension of the covered region towards higher Q 2 . A further extension towards 
higher x was possible due to better methods of calculation of the kinematic 
variables from the hadronic final statesum The improvement is due to using 
a mixture of the electron and hadron information, and momentum balance. 
For not too high Q 2 the ZEUS data reach y values as low as y ~ 0.005, 
hence the HERA data have now a small region of overlap with data from fixed 
target experiments. Such an overlap will be important for precise gluon density 
extractions. A substantial extension to smaller Q 2 was achieved by shifting the 
interaction vertex for part of the time, by analysing events with a hard photon 
radiated in the initial state, and by improving the detector components around 
the beampipe. This allowed both HI and ZEUS to reach Q 2 values down to 
about 1 GeV 2 . 



1.1 Measurements of F2 

The F 2 data from the HERAT994 run@S together with the published data 
from fixed target experimentspQtrE] are shown in figs.|| and ||. Good agreement 
between the two HERA experiments is seen, as well as a smooth continuation 
from the fixed target data to the collider data. At fixed Q 2 F 2 rises with 
decreasing x down to the smallest values of Q 2 , although the steepness of the 
rise is clearly seen to decrease with decreasing Q 2 . Similarly at fixed x 0.1 
F 2 rises with Q 2 , the rise becoming steeper as x decreases. 

The evolution of the rise with x may be demonstrated by parameterizing 
F 2 ~ x~ x at fixed Q 2 values. The result for A as function of Q 2 is shown in 
figJJ. The parameter A clearly becomes smaller with decreasing Q 2 (although 
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Figure 2: The structure function F2 as function of x for different Q 2 values, including 
measurements from HI, ZEUS, NMC, E665, SLAC and BCDMS. The solid line is a NLO 
QCD fit made using data with Q 2 > 5 GeV 2 , from HI, BCDMS and NMC. 
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Figure 3: The structure function F2 as function of Q 2 for different x values, including 
measurements from HI, ZEUS, NMC, E665, SLAC and BCDMS. The solid line is a NLO 
QCD fit made using data with Q 2 > 5 GeV 2 , from HI, BCDMS and NMC. 
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Figure 4: The exponent A as measured by the HI collaborations using fits of the form 
F2 ~ x~* at fixed Q 2 values and for x < 0.1. A-fcso shown are Ag and A 9 as a function of 
Q 2 calculated for the MRS R2 parametrizationsEj A similar plot extracted using a rather 
different technique may be found in ref.30. 
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Figucc 5: The rescaled structure functions RfF2 versus p and log(R' F F2) versus a (sec 
tcxt)cl Only data with Q 2 > 5 GeV 2 and p > 2 are shown in the second figure. The 
rcscaling functions Rp and R' F , and a s (Q) in the definition (2), are all evaluated using the 

two loop expressions in ref.12. 
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the region in x over which the function x~ x is fitted is reduced for each Q 2 
bin, which slightly enhances the effect). 

These two rises, with decreasing x (fig-||) and increasing Q 2 (fig.||) may 
be neatly combined together: at large enough Q 2 the data in the low x region 
display the universal nanJtegge rise from a flat boundary-condition predicted 
by perturbative QCDt3'til This is demonstrated in fig. |p where the rescaled 
F% is shown as function of the two scaling variables 

^^/log(^)log(^), p^log(^)/log(ggf); (2) 

the rescaling factor R' F is a simply calculable subasymptotic function, while 
Rf = e~ 2lCT R' F . The measured slope of the rise with a is 2.57 ± 0.08, to be 
compared with the QCD prediction of 27 = 2.5 fprJive active flavors. There 
exist other satisfactory parameterizations the data^JO but these do not predict 
the slope of the rise with x and Q 2 . 

An excellent fit to the data may be achieved using the full NLO pertur- 
bative QCD evolution equations, with a fitted boundary condition, as demon- 
strated in figs]| and ||. More details of such fits will be discussed in the next 
section. Models based on Regge phenomenology, but not including QCD evo- 
lution, generally undershoot the data even for the smallest Q 2 values shown. 

The increase in statistics means that it is possible to make significant 
differential measurements in the high Q 2 region for the first timet£Ej The data 
are beginning to become sensitive to the valence and sea quark distributions, 
but it is as yet too early for significant measurements of F3. Combining the 
data on charged current interactions of both experiments a measurement on 
the W propagator mass yields My/ — 82^5+3 GeV, consistent with direct 
measurements at the proton-antiproton colliders. The large Q 2 end of the 
da/dQ 2 spectrum for neutral current data is used to set limits on the effective 
scale A for contact interactions. Depending on the sign of the interference term 
with the Standard Model interactions, 95% C.L. limits lead to lower limits of 
A in the range 1.0-2.5 TeV. 

First preliminary results were presented by ZEUS of data taken with a 
small calorimeter placed 3m downstream aLthe interaction point to accept 
electrons deflected under a very small angleEZI The kinematical region covered 
is 0.16 < Q 2 < 0.57 GeV 2 and 3.10~ 6 < x < 3.10~ 5 . The data are shown in 
fig. |^ together with data from the muon experiment E665Q The latter cover 
the region 0.0DO8 < x < 0.06 and 0.2 < Q 2 < 75 GeV 2 . While the Regge 
inspired modeta predicts too low a cross section at low x for Q 2 values above 1 
GeV 2 , the agreement in this region is much better. It thus appears that there 
is a transition from the perturbative region to the Regge region at Q 2 ~ 1 



7 



ZEUS 95 BPC(preliminary) A ZEUS 94 □ E665 



0.4 
0.3 
0.2 
0.1 


0.4 
0.3 
0.2 
0. I 





L Q 2 =0.16 


L Q ! = 0.21 




Q 2 =0.27 


i 






1 Q 2 =0.34 

[ 




•• 


0^057 



1.5 
1 

0.5 




Q 2 =1.5 


L X Q 2 =3.0 

L S 4. 

— -^^^ 

i 


i- 





Figure 6: Preliminary measurement of the structure function F% as function of x for different 
Q 2 values from ZEUS, with the new beampipe calorimeter (full circles), and from the nominal 
ZEUS analysis (full triangles) and E665 (open squares). The data analysis assumes thpfe Fl = 
0. A 5% normalization error is-riot shown. Predictions for the Donnachie-LandshofO (full) 
and GR\t3 (broken) predictions are given as well. 
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Figure 7: Comparison of the results— sit different energies and triggers for the deuteron 
structure function F% from the NMCE3 The full symbols are the new (preliminary) data 
points. The inner error bars are statistical, the full ones represent the total error. 
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GeV 2 . Future HERA data will cover this region and will thus bridge the two 
presently disjoint data samples. 

New preliminary results on F2 from union-proton, and muon-deuteron scat- 
tering have been shown by the NMC collaborationEll The data taken with the 
small angle trigger in 1989 have been analysed, giving access to a lower x and 
Q 2 region than in previous NMC measurements. For incident muon beam en- 
ergies of 200 and 280 GeV, scattered muons were detected with angles down 
to 6 mrad. The proton and deuteron structure functions were measured in 
the x range from 0.0045 to 0.6 and the Q 2 range of 0.5 to 75 GeV 2 . The F 2 
for deuterons is shown in fig. [7] including the old and new datasets. The new 
NMC results compare very well with the E665 results. For the ratio of F2 
for neutrons to protons, the range covered extends to 0.001 < x < 0.8 and 
0.1 < Q 2 < 145 GeV 2 . The Q 2 dependence of this ratio shows slopes which 
are larger than expected from QCD at large x, while there is essentially no 
Q 2 dependence measured a small x. The deuteron data show indications of a 
small amount of shadowing at small x. 

1.2 Measurements of R 

New preliminary data on R — cu/o~t have been presented from the neutrino 
experiment CCFRP from NMof and from HlE3 The CCFR data cover values 
of x in the range 0.01 < x < 0.6 for the range 4 < Q 2 < 300 GeV 2 . The 
new data are in agreement with the SLAC fit through previous data. NMC 
has released new data on R measured on proton and deuterium targets, for x 
values in the range 0.002 < x < 0.12 and 1.5 < Q 2 < 20 GeV 2 . The systematic 
errors are preliminary around 0.1, flat in x, and are dominated by the 0.15% 
normalization uncertainty between the two datasets. The difference R p -R d is 
consistent with zero. 

Despite the availability of data on R for x below 0.1 from fixed target 
experiments, a substantial uncertainty still exists in the extrapolation of R 
down to the low- a; region at HERA. Furthermore at small x the quantity R is 
a sensitive probe of the gluon density in the proton. Therefore measurements 
of R at HERA are mandatory. 

While R can be measured at fixed target experiments by varying the in- 
coming beam energy (which results in a variation of y for fixed x and Q 2 in 
expression ([!]), no such opportunity exists so far at HERA. However a first 
glimpse of R at small x can be obtained by taking the measurement of F% in 
a region where the effect of R is negligible (for example y < 0.35), making 
a NLO QCD fit to these data, and then using the fit to extrapolate into the 
region of large y (for example y = 0.7). The comparison of the extrapolated 
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Figure 8: Measured longitudinal structure function Ft, by H^l for y = 0.7 as functions of 
x and Q 2 for Q 2 = 8.5, 12, 15, 20, 25, and 35 GeV 2 . The error bars are statistical errors, 
while the band is the systematic error common to all points. The full line is a perturbative 
calculation of Fjj using the quark and gluon distributions as determined by the NLO QCD 

fit to the data. 



value of i<2 with the experimentally measured cross section then gives a handle 
on R. The result is shown in fig.||. For this determination the cross section 
measurement has been extended to y — 0.7. The value of = F 2 — 2xFi is 
around 0.5 with a systematic error of about 0.2. The result is self-consistent, 
in the sense that it is compatible with the QCD prediction obtained using the 
gluon distribution found in the fit. 

Future prospects for the measuremepij of R at HERA by reducing the 
proton beam energy were also discussedC-3 A precision of about 0.3 can be 
reached on R using two beam energies (820 and 450 GeV) with the 1996 
luminosity. The measurement is again limited by systematics. 



1.3 Measurements of F£ c 

Also reported at this meeting were results on the charm contribution to F^, 
F^ ■ Hlt3 and ZEUSE3 have shown that deep inelastic scattering events with 
charm can be tagged through the detection of D° and D*^ mesons. The decays 
D° — » K~ir + and D* + — > D a n + — > K ~7r + 7r + (and charge conjugate channels) 
have been analysed. HI also demonstrated that the production mechanism for 
charm in DIS is compatible with photon-gluon fusion, the non-perturbative 
charm sea contribution being no more than a few per cent. Hence F^ c may be 
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useful as a direct probe of the gluon density in the proton. The results of the 
F^ c measurements, assuming R — 0, both for D° and D* meson production, 
are displayed in fig. ^. The errors refer to the statistics and to the experimental 
systematics. These measurements extend previous knowledge of F%° (from 
measurements by the EMC experiment E3) towards smaller values of x by two 
orders of magnitude. When combined with this earlier measurement a steep 
rise of F| c is observed with decreasing x. This is consistent with predictions 
from NLO fits. 



• H1 D* o H1 D" □ EMC 




Figure 9: as derived from the _D*+ (full dots) and the D° (open circles) cross sections 

measured by HI. Statistical (thick) and total (thin) errors are shown. The EMC data is 
also included (open squares). The data are compared with NLO predictions using GRV (full 
line), MRSH (dashed line) and MRSDO' (dash-dotted line) for m c = 1.5GeV. The shaded 
band represents the prediction from the HI NLO QCD fit to the F2 measurements. 



2 Parton Distributions and the Strong Coupling 

The success of conventional leading twist NLO perturbative QCD in describing 
all available structure function data with Q 2 ^ 1 GeV 2 , even that at very small 
x, means that we can use the data to extract input parton distributions (and 
in particular the gluon distribution) and the strong coupling constant. At this 
meeting new global fits were presented by both the CTEQ groujo (the CTEQ4 
distributions) and the MRS groupiJ (the R distributions). While the CTEQ 
starting scale Q\ generally remains at 2.6 GeV 2 , they have in the light of recent 
data also offered a distribution with Q\ = 0.7 GeV 2 ; MRS have now reduced 
the starting scale of all their latest distributions from 4 GeV 2 to 1 GeV 2 . 
Both groups now assume that at small x the input singlet quark and gluon 
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distributions behave as x~ Xs , x~ Xs > respectively, the two parameters As and X g 
being now fitted independently, rather than assumed to be equal. Note that 
these 'effective powers' must be interpreted with some care: their precise values 
are dependent on the form of the parameterization, on the choice of scheme (the 
data in fig|| are by construction in DIS scheme, while the As curve is in MS) 
and the treatment of thresholds. Nonetheless, certain qualitative conclusions 
may be drawn: both powers increase with Q 2 , As in line with that extracted 
from the HI data (see fig-||), and A g leading As- At smaller scales however, 
somewhere between 1 and 4 GeV 2 , X g goes negative and the gluon becomes 
increasingly valence-like. 

The soft (i.e. flat or valence-like) initial gluon found in these fits is of course 
perfectly consistent with the douJple asymptotic scaling seen in figj|, since this 
characterizes the non-Regge risi3 generated dynamically in perturbative QCD 
by gluon brcmsstrahlung from a flat input. It also accounts for the remarkable 
success of the GR\Eil prediction of this rise, in which valence-like inputs^! were 
evolved from a very low scale Ql — 0.34 GeV 2 : the shape of any rise generated 
dynamically from a valence-like input at a very low scale is universal (while its 
overall normalization, not calculable in perturbation theory, was fixed by GRV 
by comparison with prompt photon data, and imposition of the momentum 
sum rule)t3 

The non-Regge risA2l of the parton distributions is actually slower than 
any power of x, while faster than any power of lnl/a;, growing instead as a 
power of exp yj\alfx. In fact the data are not consistent with perturbative 
evolution from input distributions at large Qq which rise as a power, of x, since 
the growth with Q 2 of the slope of the rise is-po longer reproduced!!! (although 
in any particular Q 2 bin it is not yet easjeij to distinguish between a power 
of x and a power of lnl/a;). It is thusjujscful to consider fits with alternative 
parameterizations of the small- a; input lB 



many less inclusive 
;4 ' the analysis of the 



2.1 The Gluon 

The gluon distribution at low x is of particular interest, since it is this that 
drives the growth in the quark distribution. Althoug 
methods of determining the gluon have been proposed 
scaling violations of remains the best method at small x. Indeed the gluon 
is now fairly well determined since the F2 data are so precise, the Bjorken 
scaling violations relatively large, and driven essentially by the steep rise in 
the gluon. The results of some of the various fits are compared in fig.|i"o|. 
The HI fit, based on structure function data only (for details see ref. 4), also 
includes error bands derived from the statistics and experimental systematics. 
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Figure 10: The MS gluon density xg{x) at Q 2 = r5 GeV 2 and Q 2 = 20 GeV 2 derived from 
a NLO QCD fit of structure function data by Hip and the gluon density xg(x) »hQ 2 = 5 
GeV 2 together with results of new global fits from the CTEC£3 and MRSEj 



This error bands have been reduced with more than a factor two compared to 
the results based on earlier HERA data. Similar preliminary results became 
available from ZEUS til Note that the gluon gets steeper as Q 2 is increased: 
as explained above this is only possible in perturbative QCD if it rises more 
slowly than a power of x. 

However, care must be taken as there is no unique way of extracting the 
gluon: besides depending on the choice of scheme and the value of a g , it 
also depends to some extent on the assumed parametrization at the input 
Qg value and on the treatment of the charm thresholdE3 Near threshold it 
seems necessary to use the boson gluon fusion mechanisrp_fpr charm production 
rather than charm production via massless evolutionEj'EZI It is also difficult to 
ascertain the experimental and theoretical uncertainties. Recently the HERA 
experiments have provided a full error matrix of the systematic uncertainties 
in the data, which make a full error analysis on the gluon density possible. The 
two groups providing global fits to the data and extracting parton densities 
are strongly encouraged to take this information into account in the future. 

2.2 Calculations of F£ c , Fl, Direct Photons and Jet Rates 

An equally clean but entirely independent determination of the gluon density 
could in principle be made from an accurate measurement of—the longitudinal 
structure function F^. The results presented at this meeting^ are consistent 
with the gluon extracted from F 2 using NLO perturbative QCD (see fig. 8), 
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albeit with a large systematic error. 

A measurement of the contribution to F% which arises from charm produc- 
tion, F| c has also been advocated as a sensitive probe of the gluon density. 
It was shown recentljO that at small x the NLO QCD predictions are stable, 
with scale variations of less than ±10%, and that they offer a rather local mea- 
surement of the gluon. The dominant uncertainty in the QCD calculations 
arises from non-perturbative effects conveniently expressed by the uncertainty 
in the charm quark mass. Although the present level of precision of the data 
does not yet allow yet the extraction of the gluon from the F% c data, they 
are consistent with NLO perturbative calculations using gluons obtained from 
global fits (see fig. 9). 

The MS gluon distribution may also be determined from the 2 + 1 jet 
rate measured at HERAprta since a f-uIL-NLO perturbative calculation of the 
partonic cross-section is now available e3'Ell There is good agreement with gluon 
distributions found in global fits, the method being particularly useful in the 
intermediate x range (0.01 x Si 0.1). 

Further information on the gluon distribution, especially at larger x (x > 
0.01) can be extracted from the Tevatron pp collider data on direct photon 
and jet production. New measurements c3 of the production of single isolated 
photons show good agreement with NLO QCD calculations for photons with a 
transverse energy larger than 30 GeV. At lower Et values the data overshoots 
theory, which can possibly be explained by extra kx generated by initial state 
parton showers. The CDF inclusive single jet cross section measurement c3 
shows a remarkable excess over theory for tr an sveise ,m omenta larger than 200 
GeV. Some doubt has been cast on the assertioiOLJ that it is not possible to 
absorb the excess by adjusting parton distributions at large x: CTEQ 
found a reasonable (if not optimal) fit to both CDF and BCDMS data 
one is ready to accept a somewhat less conventional shape for the gluon at 
high x. No rise at high- Et is seen in the DO data, which however have larger 
systematic errors in this region. The precise inclusive jet data at lower Et can 
however be used to give strong constraints on the gluon in the intecaaediate x 
range, and CTEQ now include this data in some of their global fitsEa 



ve 
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2.3 The Strong Coupling 

The success of conventional NLO perturbative QCD at small x, both in ac- 
counting for the shape of F% and in consistently describing other observables 
driven by the gluon distribution, suggests that a determination of a s from 
HERA data may be possible. Furthermore, since the double scaling rise of F2 
depends strongly (and nonlinearly) on a s , but only weakly on the shape of 
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Figure 11: Contour plots of x 2 in the three orthogonal planes (Ag, A 9 ), (a s (mz), As) and 
(a s (mz), A 9 ) through the global minimum. The first eleven rnntours are at intervals of one 
unit, while those thereafter are at intervals of five unitsc3 The data are from HlQ 



the distributions at large x, precise data on F2 at small x and large Q 2 (sajz; 
Q 2 ~ 5 GeV 2 ) provide an oppjQrtunity for a relatively clean measurement liil 
The first such determinationpS using data from the 1993 run, gave a value 



a s (M z ) = 0.120 ± 0.005(exp.) ± 0.009(th.) 



(3) 



this is currently being updated using 1994 datall3 The statistical significance 
of the determination may be seen in the x 2 contour plots fig.jnj with the 
new data, ±he experimental error is considerably reduced. Despite sopfi-Lnitial 
scepticismpll the result (||) has since been confirmed in global fit£3€3 (see 
for example fig 
BCDMS resull 



It is interesting because it is somewhat larger than the 



a s (M z ) = 0.113 ± 0.003(exp.) ± 0.005(th.), 



(4) 



but is instead closer to determinations made using time-like processes, such as 
the 2 + 1 jet rate at HERAB 

The CDF inclusipi jet data with medium-S T (that is 50 £ E T ^ 200GeV) 
also seem to prefeHH the_larger value of a s , although similar data from DO 
now prefer lower valuesEjS The impact of the jet data on the value of a s in 
the CTEQ fits may be judged from fig.[l^b. MRS are now providing two sets of 
fits: Rl with a s = 0.113 and R2 with a s = 0.120. R2 provides a significantly 
better fit to both the HERA F2 data and the CDF inclusive jet data, while 
the fixed target structure function data generally prefer Rl. 
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Figure 12: \ 2 vs - a s (Mz) for global fits based on current experiments: (a) BCDMS, CCFR, 
combined HERA collider experiments, and total of DIS+DY; (b) x 2 °f CDF and DO jofcsi 
calculated using statistical errors only, using CTEQ4A and B-series parton distributions^ 

3 Higher Twists and Higher Logarithms 

Before the advent of HERA data it was widely expected! that the effects either 
of higher order logarithms of l/x (summed to give the 'hard pomeron') or 
multiparton correlations (leading to 'screening' or 'shadowing') would be very 
important at small x, and invalidate the use of conventional NLO perturbative 
QCD. At this meeting results were presented which put quantitative empirical 
constraints on the size of these effects, and explanations were put forward as 
to why they might be so small. 

3.1 Higher Twists 

The standard way of quantifying the effects of higher twists in deep inelastic 
cross sections is to rerun the perturbative NLO fits to the F 2 data using instead 

F 2 (x, Q 2 ) = F 2 NLO (x, Q 2 )(l + D 2 (x)/Q 2 ) : (5) 

an estimate of the size of the higher twist effects is then given by thejdze of 
D 2 {x) fitted in various bins in x. Fits to SLAC and BCDMS data sho"wO that 
at large x D 2 is of order 1 GeV 2 , but falls rapidly, becoming small and negative 
for x 0.3 (see fig.[l3|). Recent fits to F3 data follow the same pattern!^ while 
fits to the HERA F 2 d&tki give D 2 ~ 0.2 ± 0.2GeV 2 for 10~ 4 ^ x ^ 10~ 2 . 
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Figure 13: Higher twist corrections to F2 as deduced from BCDMS and SLAC datacj Also 
shown are the results of renormalon estimates for higher twist corrections to F2 (solid) and 

xF 3 (dashed)H 



The numbers obtained in all these estimates must be treated with caution: 
twist-4 contributions will not evolve with Q 2 in precisely the same way as the 
leading twist contribution, and leading twist evolution is difficult to compute 
reliably in the vicinity of the charm threshold. However, it seems that except 
at large x the characteristic scale for higher twist effects is no more than a few 
hundred MeV, and may indeed be much smaller. This is entirely consistent 
with theoreticaLestimates of higher twist effects in nonsinglet channels based 
on renormalonsEJ Such calculations also suggest that twist-6 contributions to 
structure functions are also very small at small x. In fact at presently ob- 
tainable energies no 'unitarization' corrections are necessary since the inelastic 
cross-section is always far below any reasonable unitarity bound. 

Of course at sufficiently small Q 2 conventional perturbation theory must 
eventually break down for all values of x, since in order to maintain finiteness 
of the cross-section Fq, ~ 0(Q 2 ), Fl ~ 0(Q 4 ) as Q 2 — > 0. However there is 
as yet no quantitative theory of inelastic lepton-proton scattering at low Q 2 . 
Fits to low Q 2 dat£3 generally have two ingredients: the x dependence given 
by Regge theory (with a pomeron intercept of 1.08) and the Q 2 dependence 
given by a combination of vector meson dominance and the photo-production 
limit. Such parameterizations give a satisfactory account of the new ZEUS 
data for Q 2 ^ 1 GeV 2 (see fig.6). The fit to the data can be maintained above 
1 GeV 2 only if the pomeron intercept is allowed to rise with Q 2 : some groups 
claim to be able to calculate this rise non-perturbatively, while others use a 
parameterization based on the rise expected from perturbative QCD. 
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3.2 BFKL Logarithms 



At small x one might expect conventional perturbation theory to eventually 
break down, even when Q 2 is large enough for a s (Q 2 ) to be small, due to 
large logarithms of 1/x. Naively one expects at n-loops terms of the form 
a"ln 2 ™ -1 1/x, but for the sufficiently inclusive quantities, such as the singlet 
component of F2, these double logarithms cancel to all orders. Single loga- 
rithms, of the form a" In™ 1/x, remain howe ver At fixed coupling they may 
be summed up using the BFKL equatioroH to give a gluon distribution which 
rises as x~ x in the Regge limit, with A = 121n2a s /7r ~ 1/2. However this be- 
haviour is far too singular to be compatible with HERA structure function 
datafJ as .explained above, a steep power-like rise would spoil double scaling. 
Recent fit£3i2j'E3 now deliver a flat or valence-like input gluon (see for example 

Many explanations have been put forward for the absence of the pertur- 
bative pomeron in inclusive quantities. In particular it is often said that since 
at present only leading order logarithms can be included, the as yet unknown 
subleading effects might somehow reduce the rate of the rise. Another ma- 
jor source of uncertainty are the contributions from the small-fcT region; the 
overall normalization of the pcrturbative rise is rather sensitive to the value 
of an infrared cutoff, and phenomcnological studies of the perturbative contri- 
bution to the total cross-section suggest that it is completely swamped by the 
soft contribution^] Nonetheless, the Lipatov approach, or some development 
thereof, remains at present our only hope of actually being able to compute 
the behaviour of structure functions in the Regge limit. Considerable effort 
has been put into numerical solutions of modified BFKL or GCFM equations, 
with some success in obtaining reasonable fits to HERA datauiaEa 

A less ambitious way of including the higher order logarithms is to ob- 
serve that at large Q 2 we can project out the leading logarithmic component 
of the BFKL equationpl retaining only an infinite (but convergent) series of 
contributions to the perturbative splitting function. Factorization, of infrared 
singularities is accomplished by the fey-factorization thcorcmpirEj a general- 
ization of the more usual mass-factorization, and likewise proven formally to 
all orders in perturbation theory. The effect of the higher logarithms is then 
interpreted not as a prediction of the form of the input distribution, but-as-a 
systematic modification of the subsequent perturbative evolution with Q 2 lp" lB 
and as such can be searched for in the data. 

There are several reasons why the resulting effects turn out to be small: 
• In the gluon channel the coefficients of the LO logarithms of 1 /x happen to 
vanish at 0(a 2 ), 0{a i s ) and O{o 5 s ), and the remaining coefficients turn out to 
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be rather small. However, although in the quark channel the splitting func- 
tions-only begin at NLO in logarithms of 1/x, their coefficients are relatively 
largeEJ Presumably the same is true at NLO in the gluon channel, although 
as yet these subleading terms have not been calculated. 

• When QL,. runs with Q 2 , the effect of the higher logarithms dies away as Q 2 
increasespl and is thus concentrated in a narrow wedge at small x and low Q 2 . 

• Much of the effect of the higher logarithms may be absorbed in the initial con- 
dition: this is the essence of factorization. When comparing calculations using 
different perturbative expansions it is thus essential to refit the boundary £on= 
dition: calculations which simply evolve from a fixed boundary conditiorEju 
inevitably over-estimate the effect of the higher-order logarithms. 

• Although the LO logarithms are scheme independent, the size of the NLO 
logarithms depends on the choice of factorization scheme. A universal (i.e. 
process-independent) reduction in the size of all NLO logarithms may be 
made by absorbing the (singular) gluon normalization factor into the input 
distributionED Furthermore, it is possible to find a scheme in which the NLO 
logarithms in the quark channel. vanished Alternatively, one can argue for the 
adoption of a 'physical' schemed in which quarks and gluons are directly iden- 
tified with F 2 and Fl at small x: the NLO logarithms are then non-zero, but 
relatively small. 

• Subleading logarithms may be chosen ixl various ad hoc ways in order to 
impose momentum conservation by handEJ In fact it is possible to choose a 
particular class of factorization schemes in which momentum is conserved at 
NLO, fixing completely theJNLO logarithms in the gluon channel in terms of 
those in the quark channclEJ Momentum conservation generally reduces the 
subasymptotic effect of the logarithms, since the subleading terms subtract 
from the leading. 

• Retaining the higher order logarithms in the large- a; region is clearly meaning- 
less: splitting functions at more than two loops will only be well approximated 
by their leading In 1 jx approximation when x is small. It is thus important 
to only include the higher loop contributions when x < xq, with rr- n a su itable 
(but essentially unknown) parameter. In practice this may be doncjO by re- 
placing In i by 9(xq — x) In 24 in all terms beyond two loops, thereby ensuring 
continuity at x — xq. As xq — > the effect of the higher order logarithms is 
then reduced systematically to zero: the empirically determined value of xq 
thus gives a useful parameterization of the constraints put on the size of the 
logarithms by the HERA data. 

The first calculations including higher-ordepJpgarithms in which all these 
considerations were properly taken into accouncaea found that with the 1993 
HERA data excellent fits could be achieved at NLO with xq = 0.1 in all 
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Figure 14: The \ 2 1^1 the fits to the HI 1994 data as a function of xq for a variety of 
factorization schemes^ Those in the left-hand figure are all MS schemes, while those in 
the right-hand figure are DIS schemes. The dotted curves are two loop calculations with 
LO BFKL logarithms added, while the dashed and solid lines are complete NLO order 
calculations in the schemes of refs.62,66 respectively. The boundary condition is refitted at 
each value of xq, and all the NLO schemes conserve momentum. 

schemes, while values of Xq as large as unity were acceptable only when LO 
logarithms are summed, or in NLO schemes where the NLO logarithms are 
suppressed. However it was also notecd that there were large scheme depen- 
dent fluctuations at small x and low Q 2 . Since the 1994 data, besides being 
much more precise than the 1993 data, also explore this region, they are much 
more sensitive to the effects of the higher-order logarithms. It turns out that 
the new data effectively exclude such logarithms: in each scheme the best fit 
to the data is obtained as xq — ► 0, and in the complete NLO calculations a 
reasonable fit is only obtained when xq is unnaturally small (see figjl^). 

It should also be noted that the NLO logarithms in the quark channel 
would lead to a suppression of the gluon distribution extracted from F 2 at 
small a; by a factor of around two, and thus to a corresponc^^uxmression of 
both Fl and F% c as compared to the two loop expectationEirEj : 0'E3 No such 
suppression is seen in the data, which if anything lie a little above the two loop 
predictions (see figs.|| and ^, also ref.72). 

In our opinion a completely satisfactory explanation for the apparent ab- 
sence of these contributions to perturbative evolution at small x has not yet 
been put forward. 

3.3 Other Logarithms 

The cancellation of the double logarithms in the singlet component of F 2 , Fl 
and F 2 ° does not occur for either nonsinglet components of F 2 or ind ee d fo r 
nonsinglet or singlet contributions to the polarized structure function ^Ljotj 
When the double logarithms are naively incorporated into, the .evolution equa- 
tions in the same way as the single logarithms were abovel, 75 H 6 1 large effects at 
small x can be generated dynamically from flat or valence-like inputs: for ex- 
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ample grows so rapidly that its first moment apparently diverges. However 
it would be premature to take these predictions at face value. As yet there is 
no fey-factorization theorem when the double logarithms are uncancelled, and 
it is not even clear whether it might eventually be possible to prove one. A 
better understanding of double logarithms will probably require first a true 
understanding of single logarithms. 

In this connection, it should also be noted that in many less inclusive 
quantities, and in particular observables related to jets, the double logarithms 
are also uncancelled^ 

All the logarithms discussed so far are ultraviolet logarithms, in the sense 
that they only become important at high energies. At lower energies infrared 
logarithms can be more important. This applies particularly to the evolution 
of structure functions at large x, where it becomes necessary to consider terms 
of the form a™ In (1 — a;). Suck.terms might be important in the determination 
of a s from fixed target datapil in the interpretation of high-i/r inclusive jet 
rate as_measured by CDF (through the comparison with evolved fixed target 
data)pl and in the evolution_qf the pomeron structure function (since the gluon 
distribution in the pomerorO seems to be very strongly peaked at large /3). 

Consider the theoretical error quoted in the determinations @ and (Q) of 
a s . In both cases the dominant contribution to the theoretical errors comes 
from a crude estimate of NNLO corrections based on a variation of the renor- 
malization and factorization scales in the range 0.25Q 2 < /i 2 < AQ 2 . The error 
in the small x determination (||) also includes an error from the effect of the 
logarithms of 1/x (though we now know empirically that the effect of such 
logarithms is very small) : the corresponding determination of the error in the 
large x result (Q) due to infrared logarithms has yet to be made. 

There are several indications that the effects of infrared logarithms on 
the analysis of fixed target data at large x may be much more greater than is 
usually supposed. Firstly, when the W 2 cut on the BCDMS_data is raised from 
10 GeV 2 to 30 GeV 2 , the value of a s increases substantially^ Secondly, for first 
moment sum rules such as the Gross-Llewellyn- Smith sum rule, which obtain 
their dominant contribution from the large x region, the NNLO correction is 
known to be_as large as 35% of the NLO, while the NNNLO correction is as 
large as 15%L3 All of these corrections have the same sign. Thirdly, higher twist 
corrections, which in the renormalon approach are considered to be directly 
related to the size of higher order corrections^ are strongly peaked at large x 
(see figjl|). 

A calculation of the three loop splitting function would go some way to- 
wards settling this issue: partial results in the nonsinglet channel, which may 
prove useful for error estimation, have been presented recentlylZa 
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4 Conclusions 



There has been considerable progress in both the range and precision of struc- 
ture function data over the last year: in particular the 1994 data on F2 are 
now final, there are new data at low Q 2 and x, and new measurements of R 
and F% c . These data have reinforced the main conclusion of the Paris meeting, 
namely that conventional NLO perturbative QCD works remarkably well for 
Q 2 <> lGeV 2 , matching smoothly onto Regge expectations below this value. 
The gluon distribution is now determined with an uncertainty of around 10% 
for x down to 10 . 

However there is now definite, though admittedly inconclusive, evidence 
that collider structure function data prefer a higher value of the strong coupling 
than fixed target data, more in keeping with high energy determinations from 
e + e~ machines. Moreover the tight constraints the new data place on both 
higher twists and, perhaps more surprisingly, on higher logarithms of X/x at 
small x suggest that perturbative evolution at high energies is still not yet 
completely understood. The new data are a continuing inspiration to theorists 
searching for a deeper understanding of perturbative QCD. 
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